Charged Scalar Fields in an External 
Magnetic Field: Renormalisation and Universal Diamagnetism 

Debnarayan Jana * 
Raman Research Institute, Bangalore 560 080, India 
(February 6, 2008) 



ON 



3 
Hi 

o 



> 
m 
in 
o 

o 

\o 
o> 



c 

o 
o 



The physical and mathematical mechanism behind dia- 
magnetism of N (finite) spinless bosons (relativistic or non- 
relativistic) is well known. The mathematical signature of this 
diamagnetism follows from Kato's inequality while its physi- 
cal way of understanding goes back to Van Leewen. One can 
guess that it might be true in the field theoretic case also. 
While the work on systems with a finite number of degrees of 
freedom suggests that the same result is true in a field the- 
ory, it does not by any means prove it. In the field theoretic 
context one has to develop a suitable regularisation scheme 
to renormalise the free energy. We show that charged scalar 
fields in (2+1) and (3+1) dimensions are always diamagnetic, 
even in the presence of interactions and at finite temperatures. 
This generalises earlier work on the diamagnetism of charged 
spinless bosons to the case of infinite degrees of freedom. We 
also discuss possible applications of the theory. 

PACS numbers: ll.10.Wx, 05.30. Jp 



I. INTRODUCTION 

A classical gas of charged point particles is non- 
magnetic, by Van Leeuwen's Theorem [0 . But quantum 
mechanically, free spinless Bose particles in a uniform 
magnetic field show diamagnetism [Q. One finds that 
the energy of the system in a magnetic field is higher 
than in the absence of the field. Simon j3j proved quite 
generally that the ground state of non-relativistic spinless 
bosons interacting through an arbitrary potential always 
increases in a magnetic field. He went on to extend this 
result by showing that the free energy in the presence of 
a magnetic field is always greater than the free energy in 
the absence of a magnetic field at all temperatures 0. 
An alternative proof of this result is given in All 
this work which deals with systems with a finite num- 
ber of degrees of freedom suggests that diamagnetism is 
a universal property of spinless bosons. In field theory 
(which describes systems with an infinite number of de- 
grees of freedom) charged spinless bosons are described 
by complex scalar fields. One might therefore expect 
that charged scalar fields would also show diamagnetic 
behaviour. In fact Brydges et. al. || have shown this 
diamagnetic inequality for the interacting scalar fields 
on a lattice. This result insisted on an ultraviolet cut- 
off. Here we want to extend to the continuum case with 
proper account of the renormalization. With this mo- 
tivation we study the magnetic behavior of scalar field 



theories in both (2+1) and (3+1) dimensions Q 

The paper is organised as follows. The first part 
deals with finite temperature free scalar field theory in 
the presence of an external homogeneous magnetic field. 
Here, we explicitly calculate the partition function and 
the free energy as a function of the applied magnetic field 
in (2+1) dimension. This expression is formally diver- 
gent. Using a suitable regularization scheme we compute 
the difference in the free energy (with and without the 
magnetic field) and obtain a finite answer. This differ- 
ence is also shown to be positive, thus establishing the 
diamagnetic behaviour of free charged scalar fields. We 
establish the renormalisation of the free energy through 
numerical and analytical calculations. Next we discuss 
the renormalization of the free energy in (3+1) dimen- 
sion. In this connection we also show the quantitative 
difference between the responses in odd and even spatial 
dimension. More specifically, we point out the drastic 
change of the behaviour of the free energy at strong mag- 
netic field in the zero temperature limit. 

We then move on to interacting scalar field theory in 
the second part. Here, we cannot evaluate the parti- 
tion function explicitly. Nevertheless we prove the uni- 
versal diamagnetism of scalar fields by assuming a finite 
momentum cutoff in the theory. If the theory is renor- 
malizable, then one can take this cutoff to infinity while 
maintaining finitcness of all physical quantities. In both 
cases the results obtained are exact. In fact we present a 
non-perturbative approach to prove this diamagnetic in- 
equality for this interacting case in an arbitrary spatial 
dimension. Before the conclusion, we illustrate few exam- 
ples of spinless bose systems and point out the variation 
of the susceptibility with temperature above the critical 
temperature in a physical way. 

Although this work is rederivation of a known result, 
we hope that it has a value beyond the purely pedagog- 
ical. Here in most of the paper we will talk about the 
response of the vacuum. The vacuum as a medium has 



1 At this point we would like to emphasize that it is ex- 
tremely naive and dangerous to conclude the response of the 
field theory from the corresponding N particle systems. For 
example, N fermions at T = show Pauli paramagnetism 
while its counter field theory QED vacuum is diamagnetic. 
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been a fruitful picture in trying to develop intuition about 
abelian and non-abelian gauge theories. In fact we point 
out several interesting points regarding the nature of the 
vacuum say for example QED and QCD in presence of 
an external magnetic field. 



II. FREE CASE 
A. (2+1) dimension 

In this subsection we calculate the free energy of free 
scalar fields in the presence of an external uniform mag- 
netic field. For ease of presentation we work here first 
in two spatial dimensions. The interesting physics takes 
place in the plane normal to the applied field. In the next 
subsection we will discuss the renormalization of free en- 
ergy in (3+1) dimension. 

Let $ be a complex scalar field which describes charged 
spinless Bosons. The Lagrangian density of a free charged 
scalar field in the presence of a constant homogeneous 
external magnetic field is given by 



C = (L> M $)*(D^$) - m 2 ($*$) 



where \i = 0, 1, 2, 



D ll = d M - ieA^ 



(1) 



(2) 



and m and e are the mass and charge respectively. (We 
set h = 1 and c = 1). Now, we write the complex field in 
term of two real fields <I>i and $2- 



$ _ $1+^2 $ » $1 - 



V2 ' 



V2 



(3) 



This theory has a global U(l) symmetry and therefore a 
conserved Noether charge Q, given by 



Q = J d 2 x (HiSa-SiIIa) 



(4) 



where 

Ui = d $i (5) 
The Hamiltonian density of the system is given by 

H = + n 2 ) + i(V*i) 2 + i(V$ 2 ) 2 + 

±(m 2 + e 2 A 2 )(<S>l + <f>l)-j-A, (6) 

where the current density j is given by 

J = -e ($iV$ 2 - $ 2 V$i). (7) 

We now suppose that the external magnetic field is uni- 
form in the x — y plane. 



We choose the temporal gauge ( Aq = 0). The constant 
magnetic field B is 



B d X Ay dy A x 



(8) 



where A x and A y are independent of t. 
The action of this theory is 



S = Jo J dr 



(9) 



where r is the imaginary time variable which runs from 
to /? (=l/(fcsT)), the inverse temperature. The action 
defined above is quadratic and so the partition function 
can be evaluated exactly. As is usual in finite tempera- 
ture field theory [[|, we impose periodic boundary con- 
ditions for Bosonic fields 



$(x,0) = $(x,/3). 



(10) 



Now, the partition function of this theory can be writ- 
ten as 



Z(B) = J 2?[nip[n 2 ] J X>[$i]X>[<J> 2 ] exp J 



dr d x 



m! —i- + in 2 ^ - $ 2 , n 1; n 2 ) 

OT OT 

+/i(* 2 ni - $in 2 ))] 



(ii) 



Here fi is the chemical potential associated with the 
conserved charge Q. The charge density (Q/A) has to be 
contrasted with the usual number density. The charge 
density refers here to the difference between the parti- 
cle density and the antiparticle density and hence can 
take any sign while the number density by definition is 
always positive. In that sense [i is not the usual chemical 
potential used in Grand Canonical Ensemble. We pick 
the gauge in which the vector potential A is (-By, 0), 
and expand the complex scalar field in terms of modes 
adapted to the present situation. These modes solve 
the Klein-Gordon equation in an external magnetic field. 
The eigenf unctions are labelled by one discrete (I) and 
one continuous p x quantum number and the spectrum 
depends on / only. In the gauge we choose, the modes 
are plane waves in the x direction and harmonic oscillator 
(i.e. gaussian) wavefunctions in the y direction. 

The spectrum is given by 



(2l + l)eB, 1 = 0,1,2 oo 



(12) 



The degeneracy of these states for fixed / is eAB/2n , 
where A is the area of the system. So, these modes can be 
thought of as quantized harmonic oscillators. Expanding 
the fields <&i and $2 in these modes the system reduces 
to a collection of harmonic oscillators with frequency u>i . 
By standard manipulations , we get the free energy 

as 
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F(B) = 27reAB^ 



1=0 



uji + — ln(l - exp(-/3(w ; - fi)) 



+ -ln(l-exp(-/3(w z +//)) 



(13) 



The first term in the square brackets corresponds to 
the zero point fluctuation of the vacuum and the other 
two terms are finite temperature contributions of the par- 
ticles and antiparticles respectively. Without the above 
formalism one can also calculate the vacuum energy by 
summing up the zero point energies for the modes of all 
the fields. Of course when we sum over all modes the 
energy actually diverges. One has to introduce a regu- 
larisation scheme or a cutoff so that the vacuum energy 
becomes finite. 

It is easy to notice that this zero point energy is diver- 
gent due to the summation of infinity number of modes 
(Landau levels). In conventional field theory this infi- 
nite zero point energy is always discarded; since it can 
be reabsorbed in a suitable redefinition of the zero of en- 
ergy. This is justified in the sense that the infinite zero 
point energy is unobservable. However, the change in 
zero point energy caused by the external constraint is fi- 
nite and observable. So, according to Casimir's idea, 
the physical vacuum energy can be defined as the dif- 
ference between the zero point energy corresponding to 
the vacuum configurations with constraints and the one 
corresponding to the free vacuum configurations. This 
definitions must be supplemented in general with a reg- 
ularisation prescription in order to obtain a finite final 
convergent expression. We will first compare the free en- 
ergy of the system with and without the magnetic field 
at zero temperature. 



1. Zero Temperature Case 

The free energy of the system in presence of the mag- 
netic field at zero temperature is given by 



F (B) = 2nA eB^T 



(14) 



where uf = m 2 + (2l+l)eB. Obviously, this sum diverges. 
In order to obtain a finite answer we need to impose a 
cutoff L in the sum (H). Then the free energy becomes 



F (B,L) = 2ttA eB^ uii. 

1=0 



(15) 



The free energy in the absence of the magnetic field at 
zero temperature is given by the divergent expression 

F o (0) = 2-kA f pdp y/bZ+pft + m* (16) 



We regularize this 
Then the free energy 

Fo(0,A) = 2-kA 



expression 



by imposing a cutoff A. 
(ilf) becomes 



pdp 



(pI+pI) 



(17) 



In order to compare the free energies in equations (|l5| ) 
and ([IT]), we choose the cutoffs L and A in such a way 
that both systems have the same number of modes. Such 
a procedure can be justified on physical grounds if one 
imagines that the magnetic field is turned on adiabati- 
cally j8|]. In the Appendix A and B, we show that this 
mode matching scheme is the correct regularisation pro- 
cedure. More specifically, in Appendix A we compare 
this scheme with another scheme which looks apparently 
correct and in Appendix B we justify this mode matching 
scheme from the field theoretic point of view. 

Counting the modes upto the L-th Landau level we 
find 



2ttA eB^ 1 = 2nA eB(L + 1) 



(18) 



2=0 



Similarly, for the momentum cutoff upto A we get the 
modes without the magnetic field as 



2-kA I pdp = 7r AA 2 
Jo 

Equating these gives us 

A 2 = 2 eB (L+l) 



(19) 



(20) 



Now, the free energy in absence of the magnetic field 
depends on magnetic field through the momentum cutoff 
and is given by 

rHB) 

F (0,L) = 2ttA pdp vV + m 2 (21) 

Jo 

The difference between the two free energies is given by 
AF(B, L) — F (B, L) — F o (0, L) (22) 

We define f(B) = F (B, L)/2nA, f(B) = F a (0,L)/2nA 
and Af(B) = f(B) — f(B). Numerically evaluating these 
sums one can show that for finite L, Af(B) the difference 
between two large quantities is positive. As the cutoff L 
goes to infinity, Af(B) becomes the difference between 
two infinities. In this limit we find that Af(B) tends to a 
finite value. Thus, the susceptibility at zero temperature 
in the relativistic case is non-zero. This vacuum suscep- 
tibility can be interpreted as due to virtual currents. 

We now show analytically that AF(B) is positive i.e. 
the vacuum is diamagnetic. Note that 



Af(B) = f(B)-f(B)=J2 ai(B,m) 



(23) 



1=0 
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where a;(i?,m) is given by 



ai{B,m)=eB yjm 2 + {21 + l)eB- 



da y/m 2 + 2(1 + a)eB 



(24) 



Introducing a dimcnsionless quantity p = ^| the above 
equation becomes 

f-i 

ai{p) = P 



^l + (2l + l)p - / da yjl + 2{l + a)p 
Jo 



(25) 

The positivity of ai (p) for each I can be proved geomet- 
rically. Defining Zi = (1 + 2lp)/2p and f(a) = y/zf+a, 
the coefficient ai(p) can be rewritten in terms of Ci(p) as 

ci(p) = -^f 2 = mm f Q da f(a). (26) 

Since, the function f(a) is convex, the area under the 
tangent drawn at a = 1/2 is greater than the area under 
the curve (see the figure below). 
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Figure 1: The full line is the curve f(a) = + a. The 
dashed line is the tangent to the above curve at a = .5. It 
is straightforward to see that the area under the curve is 
less than the area under the tangent. Also the area under 
the tangent is the same as the area of the rectangle. The 
area of this rectangle is given by /(1/2) x 1 = /(1/2). 
Hence, the positivity of a(.l) is proved. This can be gen- 
eralised to any positive value of z;. 



This shows that ci(p) is positive. To show the conver- 
gence of the sum (^3|) we note that 



da f(a) < 



mm 



/(o) + /(i) 



(27) 



Now, applying mean value theorem twice one can easily 
show that 



daf(a)< — 1 3/2 
o 16{zi + a) 



(28) 



Thus the coefficient ci(p) is positive for each I and the 
sum converges, hence the diamagnetic inequality is es- 
tablished. 



2. Massless limit 

In this section we want to discuss the behaviour of 
the leading term of the free energy in the massless limit 
and its consequences. It is easy to see from the zero 
temperature free energy that the magnetisation in the 
zero mass limit is given by 



M(B) ~ - VB (29) 
So, the susceptibility in this zero mass limit is given by 

X(B) \= (30) 



B 

which diverges |0] as B goes to zero. This divergence 
of the susceptibility is reminiscent of the fact that the 
free energy F(B) ~ B 3 ^ 2 . This variation of the free en- 
ergy in the massless limit can also be understood from 
dimensional arguments as follows. Since we are working 
in natural units Ti = 1 and c = 1, then [m] ~ 
Then the free energy density (i.e. per unit area ) varies 
as [i] -3 - However, the dimension of B is [£]~ 2 - So, the 
massless limit restricts the free energy density variation 
with magnetic field B to B 3 / 2 only. This feature of the 
susceptibility has already been noticed in the magnetised 
pair Bose gas |1C[| . This divergence of the susceptibility 
can be compared with that of ideal charged bose gas at 
the condensation point Q . The magnetic susceptibility 
is defined by 



M = xB 



(31) 



where B is the applied field; what is actually computed, 
however, is a quantity \' defined by 



M = X 'B' 



(32) 



where B' is the "acting" field. The relation between B 
and B' requires special attention; in fact if the "acting" 
field is identified with the average microscopic field then 
we have for this case 
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B' = B + 2ttM 



(33) 



This equation immediately gives a relation between \ 
and x' as 



X' 



1 - 2ttx' 



(34) 



This shows that on approaching B — ► limit, where 
X 1 — oo, — o - so that the permeability tends to 

zero. So, in this case the external field will be totally ex- 
pelled. This happens because of large number of virtual 
particle and antiparticle produced in the ground state 
so that the overall diamagnetism of the system is high 
enough to totally expel the external field. That the mag- 
netisation is non-analytic in the limit i? — > is the actual 
signature of a phase transition as a function of external 
field strength. Note that this behaviour is equivalent to 
the behaviour of the free energy in strong magnetic field 
limit. This can be understood the way to and B appears 
in lo — ^ which is the only energy scale in the problem. 
Hence, one should expect the same behaviour as to — * 
or B — ► oo because in both situations higher Landau lev- 
els become less and less occupied. In fact in this situation 
one can construct the theory in the lowest Landau level. 

Another interesting point is that there exists a crit- 
ical magnetic field below which magnetic field will be 
totally expelled. This critical field can be estimated as 
follows. The effective magnetic field can be defined as 
B e ff = B + 2ir M . Here, the magnetisation M varies as 
— const\f~B. Therefore, there exists a critical magnetic 
field B c where the effective field B e ff vanishes. 



3. Finite Temperature Case 

Now, for the finite temperature case one can regulate 
the free energy through the same mode matching regu- 
larisation method. Finally, one can write down the free 
energy difference in dimensionless form as before 



AF(B) = F(B) - F(Q) = ]T k{p, 8, 



where, 



1=0 



g(p,l, 1/2) - / da g(p,l,a) 
o 



(35) 



(36) 



The dimensionless variables are defined as 8 — (3m and 
The coefficient g(p,l, a) is given by 

g(p, I, a) = log (l - exp(-S(y/l + 2{l + a )p - 0)) + 

log (l - exp(-5(Vl + 2(1 + a)p + 0)) • (37) 



Now, defining z; 



S 2 (l+2lp) 
2p 



we can write the equation 



g(p, I, a) = log (1 - exp(-(v / z/ + a - 0)) + 

log (1 - exp(-(V5FFS + 0)) . (38) 

The function g(p, I, a) is convex, so the zero temperature 
argument applies unchanged. It follows that the free en- 
ergy satisfies the following inequality 



F(B) > F{0) 



(39) 



Thus the response of the system to the magnetic field 
will be diamagnetic. 



B. (3+1) dimension 

In this subsection we would like to renormalize the free 
energy in presence of an external uniform magnetic field 
in (3+1) dimension. For the sake of repetition we do not 
present here the derivation of the free energy in (3+1) 
dimension which can be easily be done in the framework 
presented earlier for (2+1) dimension case. In fact the 
final result can be understood easily from the dispersion 
relation (energy spectrum) for (3+1) case which is given 
by 



Vp^ + rn 2 + {21 + TjeB, 1 = 0,. ...oo (40) 



Thus, we get the free energy of charged scalar fields 
in 3d in presence of a constant homogeneous external 
magnetic field as 



°° />oo 

F{B) = 27reVB^ / dp. 



iln(l- e -^ l+ ^) 
P 



(41) 



From the above equation (^lj) we notice that the zero 
point free energy is given by 



°o POO 

F(B,m) = 2neBV^2 dp z 

i=o 



a/to 2 +p 2 z + (2l + l)eB 



(42) 

We regularise this free energy by the prescription given 
in jl2| which was used in case of QED in an external 
magnetic field. In principle, one can use the prescription 
through a novel mode matching already developed for 
(2+1) dimension to regularise this free energy . But we 
use this regularisation scheme to compute quantitative 
evaluation of the regularised free energy in two interest- 
ing limits in the zero temperature case. As it is clear 
from (2+1) dimension case that we cannot evalute the 
response function quantatively that's why we choose this 
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scheme for that reason. 

In field theory there are various regularisation schemes 
for example Schwinger's proper time regularisation 
scheme fj"3| or invariant regularisation scheme due to 
Pauli and Villars [[L4J. Through this simple regularisa- 
tion scheme we demonstrate the renormalisation of field 
strength and charge as shown in the context of vacuum 
polarisation 



1. Zero Temperature Case 

The divergence of the above integral (Eq. (|42|)) is elim- 
inated by subtracting the value of the sum when B = 0. 
Hence, to carry out this renormalisation procedure it is 
appropriate to calculate first the convergent expression. 
Defining the free energy density / — F/V we can write 



df 
d(m 2 ) 



2-neB 



E 



dp z 



yjm 2 +pl + (2l + l)eB 



(43) 



d 2 f 

Differentiating once again we get = g^iy 



= -ireB 



Y^(m 2 +p 2 z + (2l + l)eB)-^ 



.1=0 



dp z 



(44) 



Notice that the integral is now absolutely convergent. 
Performing the integral and the sum, we obtain 







-ireB 



E 



i 



(21 + l)eB 



ireB I i 
>o 



-~(m 2 +eB)r] 



1 



l_ e -2eB V 

which reduces to a simple integral form given by 
neB f°° dr) e-" 1 ^ 



dri, (45) 



= 



ireB 



smh(eBr)) 
dr/e-™ 2 " g{eBrj), 



(46) 



where g(eBr/) — 1/ sinh(e.B?7). Now to get back the free 
energy we have to integrate twice. Integrating once we 
get with a undetermined constant 



0i = j d{m 2 ) 
ireB 



(-1/rj) e- m » g(eB V ) dn + C 2 . (47) 



Similarly integrating once more one gets back 



02 = 



0i d{m 2 ) 



neB 



1 



to T 
+ C 2 m 2 + Ci 



e~ m v g(eBr)) dr) 



(48) 

These two undetermined constants depend only on B but 
not on to 2 . Now it is easy to take B — > limit in the 
above formula. This gives 

02(5 = 0) = ! / -e~ m2rid V (49) 
1 Jo V 

Here we have assumed that two constants C\ and C 2 
vanish as B — > 0. Therefore the difference between the 
two free energies can be written in an integral form as 



A/ = - 
J 2 



1 

if 



dr) 



eB-q 



smh.(eBrj) 



- 1 



+ C 2 to 2 + Ci 



(50) 



Till now the constants C\ and C 2 are undetermined. 
These are determined from the following considerations. 
Parity and dimensional analysis give us a particular form 
H of A/ given by 

' B 2 ' 



A/ 



4 

to g 



(51) 



From equation (|5l]) it is clear immediately that there 
will be no odd terms of to 2 in the expression for A/. 
Hence, C 2 is zero. Now to find out the coefficient C\ 
we follow the following prescription. Notice that the ex- 
pansion of A/ in B 2 begins with a term linear in B 4 . 
Therefore, a term in B 2 would alter the coefficient in the 
Lagrangian Lq = This essentially modifies the defi- 
nition of charge. Hence, the elimination of B 2 term thus 
corresponds to a renormalization of charge. 

Expanding sinh(e_B?y) in leading order upto B 2 we get 
e 2 B 2 rr 2 e~ m ^ drj 



Ci = - 



Rescaling m 2 ?? 



o W 
rj and defining b = 



dr) 



br) 



(52) 

-^1 finally we obtain 

777. J 

(53) 



1 



1 7.2 2 
- 7) 

6 ' 



r)° [sinh(6?7) 

This renormalized energy difference A/ is one of the 
main results of the paper. Notice that the difference of 
free energy A/ goes to zero as B — > as it should be 
since we have subtracted the zero field free energy. One 
can easily show through numerical integration that this 
free energy is positive definite and it converges. Before 
we end this subsection we want to make a remark on this 
regularisation scheme. A careful reader will immediately 
notice that this scheme will not work in (2+1) dimension. 
This shows the important role played by the p z integral. 
In the next two sections we consider two interesting lim- 
its of the free energy difference. 
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2. Small Magnetic Field Limit 

At small magnetic field {brj being very small) we can 
approximate the free energy difference upto the leading 
factor in b 4 as 



A/ 



7irm b 

720 , 
77re 4 B 4 
720m 4 



e q r\ drj 



(54) 



Thus the free energy difference at small magnetic field 
limit goes as B A . This is obvious from the fact that 
we have subtracted the leading divergence of B 2 term 
from the free energy. So, the next order term in the free 
energy difference is B A according to symmetry. Notice 
that the free energy difference is positive definite and 
thus establishes the diamagnetism of the charged scalar 
fields at zero temperature. 



3. Strong Magnetic Field Limit 

In this section we would like to address various re- 
sponses of the system in the strong magnetic field limit. 

In a strong magnetic field (b 3> 1) rescaling bjn—t r\, 
we get the free energy difference from equation (|53|) 



A/ 



-n/b 

rj 3 



V 



sinh(ry) 



1 

i + - 6 v 



(55) 



When b 3> I, the important range in this integral is 1 <C 
rj <C b, in which e~ v / b w 1, then terminating the range 
of integration (with logarithmic accuracy) at rj w 1 and 
rj w 6, we get 



4^2 



A/ 



mn b 
12 



ln& 



(56) 



It is important to note that in the massless limit the free 
energy in 3d turns out as 



J 12 \ eB 



(57) 



Bo 



HL ( 3 



(59) 



where a 



is the fine structure constant. Before we go 



to discuss the response of the system in an external mag- 
netic field we would like to make some comments on the 
response of electric field alone. Since an external electric 
field always gives a non- vanishing probability for pair cre- 
ation, it is more convenient to discuss the response to an 
external magnetic field and calculate the magnetic per- 
meability \x rather than e. Then we can deduce e from 
fi if we assume that the vacuum is Lorentz- invariant. If 
that is the case then we should have 



e/i = 1 



(60) 



Depending on the value of /i < 1 or // > 1 one can classify 
the nature of the vacuum. 

From equation ( p6| ) it is easy to calculate the magneti- 
sation and the susceptibility in this limit. The magneti- 
sation is given by 



M(B) 



TreB m 2 ire 2 B 



12 



6 



In 



eB 

2 



(61) 



Another interesting point before we go to estimate the 
susceptibility is that there exists a critical magnetic field 
below which magnetic field will be totally expelled. This 
critical field can be estimated as follows. The effective 
magnetic field can be defined as B e ff — B+4ir M. Here, 
the magnetisation M varies according to equation (|6l]). 
Therefore, there exists a critical magnetic field B c where 
the effective field B e // vanishes. This critical magnetic 
field B c is given by the following equation 



In 



eB, 



6 



nem 
12 



g{e,m 2 



and in a more simple form 



B c = — exp(g). 

e 



Therefore, the susceptibility is 



(62) 



(63) 



where A is some UV cutoff in the theory. The form im- 
mediately suggests the equivalence of massless limit and 
large magnetic field limit. This is quite obvious from the 
way B and m occur in u> = ^— ■ which is the only energy 
scale in this problem. If one compares with Lq then one 
finds 



L 



2n 2 e 2 , [eB 
3 V m z 



(58) 



which implies that the radiative correction to the field 
equations may become of order of unity in strong mag- 
netic fields given by 



X(B) 



6 



hi 



(64) 



Notice that in 3d, the zero field susceptibility has a pe- 
culiar logarithmic form in the free energy. Even for 
QED [|l7j one also gets same logarithmic behaviour of 
the free energy (and hence in the susceptibility) in this 
strong magnetic field limit. For the sake of comparison 
we notice that 



X3d(B)~-ln(B) 



(65) 
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This interesting behaviour of the susceptibility in the 
massless limit is reminiscent of the fact that the free en- 
ergy variation in two cases as 



hd{B) 

hd{B) 



B 3/2 

B 2 ln(S) 



(66) 



Also it is interesting to note that in both cases the free 
energy is a non-analytic function of B. 
One can also calculate the critical magnetic field for 
which the effective permeability becomes zero. The per- 
meability is given by 



H{B) = 1 + Attx(B) 



, 2 27T 2 e 2 1 f eB 

1 — 7r em — In — ^ 

3 \ ra z 



Thus, the critical magnetic field B\ is simply 

™2 



B, 



■ exp(d(e, m )), 



(67) 



(68) 



where the function d(e, m 2 ) is determined from the equa- 
tion 



d(e, 771 ) 



27r 2 e ; 



(l — 7rem 2 ) 



(69) 



Thus we have three critical values of magnetic field for 
which (i) B e ff vanishes, (ii) p{B) vanishes and (hi) log- 
arithmic term in the susceptibility goes away and x{B) 

2 

becomes — 1ze ™ . One thing to notice here that all the 

2 

critical magnetic fields are of order — and these crit- 
ical fields should be compared with the field on which 
the classical theory of electrodynamics will break down. 
This classical magnetic field can be estimated from the 
comparision of two length scale 



e 
rn 



m 
eB 



This gives the critical field 



On comparison we notice that 

B r ] ri a. 



1 



B r 



(70) 



(71) 



(72) 



This is true for all the other critical fields also. Thus, the 
ratio between the classical field and the critical field is of 
the order of fine structure constant. 

Next we would like to compute the dielectric constant 
e(r) from the effective permeability. Since we are working 
in natural units K = 1 and c = 1 then the dielectric 
constant e and the permeability p is related by ep — 
1. However, one can show from the effective action 



formalism that under certain circumstances one can still 
use the relation in presence of an external magnetic field 
though the physical situation is not Lorentz-invariant. 
Hence, the effective permeability p(B) should be related 
to the dielectric constant e(r) by the relation given by 



e(r) 



1 



KB) 



leB->l/r- 2 



This immediately suggests that 
e(r) 



1 - 7r 2 em 2 - 



W— ) 



(73) 



(74) 



From equation ( p7| ) we notice that p(B) < 1 which means 
that e(r) > 1. This is easily seen from the above equa- 
tion. Hence, the charged scalar field vacuum (diamag- 
netic) is screening one. Before we go to finite tempera- 
ture case scheme we would like to make a comment on the 
vacuum energy of QCD in an external magnetic field. In 
fact through this regularisation scheme one can explicitly 
show the paramagnetic nature of this spin-1 field |C3|. 



4- Finite Temperature Case 

Now, for the finite temperature case one can regulate 
the free energy through the same mode matching regu- 
larisation method used in (2+1) dimension. Finally, one 
can write down the free energy difference in dimensionless 
form as before 



AF(B) = F(B)-F(Q)=J2 <MM,C,p) 



(75) 



where, 



di(p,S,() = 



dp z 



g(p,l, 1/2)- I da g(p,l,a) 
o 



(76) 



The dimensionless variables are defined as 6 = (3m and 
P = 0fi. 

The coefficient g(p,l,a) is given by 
g(p, I, a) = log (l - exp(-S{y/l + 2(l + a)p - C))) + 

log (l - exp(-<S(v/l + 2(1 + a)p + C))) . (77) 

we can write the equation 



Now, defining zi 



(l+2ip) 
2p 



g(p, I, a) = log (1 - exp(-(V^ + a - 0)) + 

log(l-exp(-(V^ + ^+0))- (78) 

The function g(p, I, a) is convex, so the zero temperature 
argument applies unchanged. It follows that the free en- 
ergy satisfies the following inequality 
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F(B) > F(0) 



(79) 



This proves the diamagnetism for charged scalar fields at 
finite temperature. It would be interesting to look into 
the various limits of this free energy difference. But here 
we do not address these limits at all. 



III. INTERACTING CASE 

In this section we want to extend the diamagnetic in- 
equality to the self-interacting field theory case includ- 
ing the dynamical interaction between scalar fields. We 
present here the proof of diamagnetic inequality for a 
generalised d dimensional case. The partition function 
of this charged self-interacting field theory in the pres- 
ence of the magnetic field can be written as 



Z(B) = J f V[<f>]V[<S>*]exp(-S '($,$*, A)), 



(80) 



where the action S is defined as 

S = J J d d x dr [(£>,,$) (L>^$)* + m 2 ($*$) + V(**$)] • 

(81) 

The action is not quadratic and Z(B) cannot be eval- 
uated in closed form. Nevertheless, we can show that the 
response of the system to an external magnetic field is 
diamagnetic. Since the formal expression for the parti- 
tion function may not exist (the integrals may not ex- 
ist) we impose a cut off in momentum space. The func- 
tional integral in (80) signifies that one only integrates 
over those field configurations whose Fourier transforms 
have support within a sphere of radius Ao in momentum 
space. The partition function then explicitly depends on 
Ao. We do not explicitly indicate the A and \i depen- 
dence of Z(B, Ao,/z) below. 

We divide the action into two parts Sq and Sj nt , where 
So is the action in the absence of the external field. 



S — Sq + Sj nt , 



(82) 



where 



s °"/ / ddxdT [(^ $ )(< 9Ai *)*)+ m2 ( $ * $ ) 

+V^(<r$)] , (83) 



Si„t = j J d d x dr [-ie(d^)(A^*)+ie{A^)(d^*) 

+e 2 (A- A)($$*)] . (84) 

Notice that exp(— So) is a positive measure on the 
space of field configurations. The ratio Z(B)/Z(0) 



can therefore be regarded as the expectation value of 
exp(— Sint)- Since exp(— S; n t) is an oscillatory function 
whose modulus is less than or equal to 1, we conclude 
that 



Z(B) 
Z(0) " 

This implies that 



=< exp-S mt >< 1 



F(B) > F{0) 



(85) 



(86) 



This result is an exact and non-perturbative one. 
Hence, it is more general and strong compared to per- 
turbative one. 

In this derivation, we have not assumed any form for 
the vector potential. So, the result derived above is true 
for both homogeneous or inhomogeneous magnetic fields 
of any strength. Since (3 is arbitrary, the result holds 
at all temperatures. The argument presented here works 
for any arbitrary interaction V(Q*Q) ( Generally, it is 
assumed that y (<&*$) is a smooth function, for instance, 
a polynomial). 

Upto now we have considered the cases of charged 
scalar fields interacting through a potential. It is also 
possible to consider interaction mediated by a dynamical 
electromagnetic field A^. The fields in the system are 
now $ (charged scalar fields) and A^. If one applies an 
external magnetic field A ext then the full Lagrangian is 
given by 



1 



F 2 + (£>,,$)* (£>"$) - rn 2 (<I>*$) - !/($*, $) 



(87) 



where D,, 



ieAf* - 



djiAi, 



'is "/I • 1 >,, and l) tI , — (V^.l,, 

The argument given above can be modified as follows. 
The definition of So changes slightly while Si n t remains 
the same. 



So = 



and 



d d x dr 



-If 2 



2 ($*$) + !/($*$)] , 



(88) 



S int = / J d d xdT[-\e(d^){A^*) 

+ie(A Al $)(cK 1 $*) + e 2 (A • A) ($$*)] . (89) 

Again one can repeat the same argument to establish 
the diamagnetic inequality by noting that exp(— So) is 
a positive measure and the ratio Z(B)/Z(0) as an ex- 
pectation value of exp(— Si n t). This universal inequality 
follows from basic principles and does not depend on the 
details of the interaction. 
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A. Comment on the regularisation scheme in 
interacting case 

The effective potential V / (<E>*<E>) is gauge invariant. 
What we have done is that we have suitably redefined 
the measure of the integral through the inclusion of 
self-interaction. Thus now, the bare measure has been 
changed to so called "dressed" measure. This has the 
advantage that the effects of self-interaction and mutual 
interaction have been separated. Also In field theory one 
would also require that the interaction P r ($*$) be renor- 
malizable. In (d + 1) dimensions this would restrict the 
interaction (<!>*<I>)™ 1 only. From simple power counting, 
one can notice easily that the required value of m is j^j . 
In condensed matter physics, where a natural cutoff ex- 
ists, higher order powers of (<&*$) may also be present. 

Under normal conditions kgT <C m we expect the self 
interaction to merely affect the particle's mass. This is 
the same concept one uses in quasiparticle picture of Lan- 
dau Fermi Liquid theory |^(| . The effect of all the inter- 
actions can be put in the definition of the effective mass of 
the quasi particles. Therefore, the error committed would 
be small if we were to ignore the self-interactions while 
keeping each mass at a fixed observalble value. Hence, 
we drop the superscript "dress" and work with the ordi- 
nary Wiener measure. 

Notice that the bare mass enters the Wiener measure. 
However, in the expression of the Sj n t neither refers to 
nor uses the bare mass; it never enters the interaction of 
the scalar particle with the electromagnetic field. On the 
othcrhand, it is important to realize that the bare mass m 
differs from the physical or "effective" mass m e ff. The 
differnece is induced by the interaction. Adjusting the 
bare mass so as to gurantee that m e // agrees with the 
observed mass of the particle constitutes the central issue 
of the renormalization program. It is generally believed 
that this issue ought to be resolved outside the perturba- 
tion theory. In this paper we do not attempt to answer 
these questions. 

In the interacting case it turns out that the partition 
function and hence the free energy cannot be evaluated 
in closed form but still one can show the diamagnetism 
of charged scalar fields. This shows that diamgnetism of 
charged scalar fields is universal i.e. robust in the sense 
that it is true for all dimensions and independent of the 
interaction between the scalar fields. One more comment 
regarding the meaning of universal behaviour of this sys- 
tem towards an external magnetic field. The proof of 
the diamagnetic inequality in the interacting case does 
not assume the form of the vector potential. Thus, in 
that sense the theory is true whether the applied mag- 
netic field is homogeneous or inhomogeneous. The word 
universality also reflects this fact that the behaviour of 
the system towards the magnetic field is independent of 



the nature of an applied magnetic field. We would like 
to emphasize one point regarding our theory which is as 
follows. If a phase transition takes place changing the 
drastic nature of the system, then of course our theory 
fails. As long as a phase transition does not take place 
that alters the bosonic behaviour, this theory can predict 
the response of the system in an external magnetic field 
as diamagnetic. 

IV. APPLICATIONS 

It is evident from the previous sections that though the 
results obtained are quite general and important from 
the theoretical point of view but fail to connect directly 
with the experiment. What we have shown is basically 
the nature of the response of a system under an external 
magnetic field. Thus, though the results are of interest- 
ing one but fails to compute explicitly the susceptibility 
as a function of magnetic field, mass and temperature. 
Only in zero temperature case we have been able to com- 
pute the susceptibility in two limits. However, one can 
resort to perturbation method to compute approximately 
the response of the system in case of finite temperature. 
In this section we would like to discuss the possible sit- 
uations which could be well-connected with the theory 
presented. From an experimental point of view one is 
more interested in the quantitative variation of the sus- 
ceptibility with the temperature. Therefore, in this sec- 
tion we will mainly address the temperature variation 
of the diamagnetic susceptibility assuming an effective 
Landau-Ginzburg free energy functional apart from giv- 
ing some illustration of spinless bose systems. Notice 
that the order parameter in Landau-Ginzburg model is 
a scalar one which has a characteristic U{\) symmetry. 
Only the mass parameter in Landau-Ginzburg model has 
a different meaning in contrast to the theory presented 
here. 

One obvious example spinless bose systems in the lab- 
oratory is the Cooper pair formation in superconductors 
which shows perfect diamagnetism (known as the Meiss- 
ner effect) below the critical temperature. Cooper 
pairs also exist in Neutron Stars where the magnetic 
field is very high compared to any laboratory field. Of 
course, the operators which create and destroy Cooper 
pairs are not strictly Bose operators, so this is only an 
analogy. In fact one can also calculate the diamagnetic 
susceptibility of the Cooper pairs above the critical tem- 
perature, more specifically the temperature variation of 
the susceptibility. The diamagnetic susceptibility above 
the critical temperature has been observed recently in 
High-T c compounds [^3). In literature there exists a 
model |24j] which does predict that above the critical 
temperature all or at any rate some of the carriers are 
bosons. Here, of course, as has been pointed out 
that the definition of the "critical temperature" bears a 
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different meaning in contrast to usual one. The "critical 
temperature" in that context refers to the temperature at 
which the gas of bosons become wholely non- degenerate. 
However, in the following we do not address all these is- 
sues regarding the mechanism of formation of bosons. 

With the usual Landau- Ginzburg effective field the- 
ory [p6| it turns out that the diamagnetic susceptibility 

4 — d 

in d dimension Xd ~ —(T — T c ) ~ . Without going 
much details into this calculation (^tJ, the above pecu- 
liar variation of the susceptibility with temperature can 
be explained as follows. Following the argument given for 
3d p8| , we notice that above T c , droplet of Cooper pairs 
will grow and decay as a result of thermodynamic fluc- 
tuations. Their mean radius is approximately equal to 
£ which is the coherence length of the Cooper pairs and 
phenomenologically the simplest variation of this coher- 
ence length is taken as £ 2 ~ (T — T c ) . The "mass pa- 
rameter" in this theory is inversely proportional to this 
coherence length. Therefore, the amount of energy re- 
quired to produce a droplet is given by 



SE ~ i x £ d x |0| 2 



(90) 



Here \4>\ 2 is the density of the Cooper pairs in the frame- 
work of phenomenological Landau-Ginzburg free energy 
functional. This energy must be equal to the thermal 
energy fcsT c and hence 



ksT c 



(91) 



Now consider the expression for the diamagnetic suscep- 
tibility of atoms in scaled form as 



Xd 



-(T - T c 



(92) 



In otherwords for a general variation of the coherence 
length £ 2 ~ (T — T c )~ v the susceptibility variation turns 
out as Xd ~ — {T — T c )~ u ( A ~ d ^ 2 . However, as it is ev- 
ident for dimension d > 4 the above formula does not 
make any sense; but this formula of course correctly re- 
produces the expected temperature variations for 2d and 
3d. In fact, in 3d the exponent of the susceptibility ( 
1/2 ) has been observed in the experiment |23|. Notice 
that inspite of the strong correlation among the Cooper 
pairs in the system the pairs behave as if they are non- 
interacting. Thus, within this simple Landau-Ginzburg 
approach one can give a strong hint towards the role of 
dimensionality of the system in High-T c material. 

It has been shown explicitly by J. Daicic et. al. |h]] 
that the magnetised pair Bose systems are relativistic 
superconductors. These systems are not covered by pre- 
vious analysis [[l 29 |^] which apply to non-relativistic 
quantum mechanical systems. Pions would be suitable 
candidates for application of this theory with ir + and ir~ 



regarded as the particles and antiparticles. The choice is 
also motivated by the fact that these pions are massive 
(rac 2 = 139.5673 Mev), obey Bose-Einstein statistics and 
that they possess no spin. It is also well known that hard 
core bosons pp| in any dimension on any lattice show a 
preference for zero flux. 



V. CONCLUSIONS AND PERSPECTIVES 

The response of a system to an electric field is com- 
pletely different from its response to a magnetic field. 
The basic difference between the responses of a system 
on application of an electric field or a magnetic field lies 
in the Hamiltonian of the system. 

The Lagrangian of a system in the presence of an elec- 
tric field can be written as 



C = (D $)*(D <Z>) - (V$)*(V$) - m 2 ($*$) 



where 



D = do - ieA 



!/($*$) 
(93) 

(94) 



For statistical mechanics to make sense, the Hamiltonian 
H must be independent of time. Then it follows that 

n = (n*)(n) + (v$)*(v$) + ™ 2 ($*$) + 

-ie[(ir)(4,$)-(II)t4o**)] (95) 



The electric field appears in the Hamiltonian through 
the linear vector potential Aq term. Now, from finite 
temperature second order perturbation |3l],|3^] theory, 
one can show easily that the free energy of the system 
always decreases with the electric field. Hence, the di- 
electric susceptibility is always positive in thermal equi- 
librium. 

But in the case of a magnetic field the Hamiltonian 
contains both linear and quadratic terms in A. The net 
effect of an applied magnetic field is not a priori clear. 
However, as our analysis makes clear, for charged scalar 
field theories the net effect is always diamagnetic. 

In case of spinless bosons, there is no Zeeman term 
coupled with a magnetic field and hence the system con- 
sisting of spinless bosons always has higher energy in a 
magnetic field than without the magnetic field. It has 
been already pointed (|^,^,|| out that there is no cor- 
responding theorem for fermions. So, in case of fcrmions 
having spin the general tendency is to show paramag- 
netism Q. 

Before we end we would like to comment on a re- 
cent work in the literature fl34fl on ultrarclativistic hot 
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scalar plasma. The results of their perturbative treat- 
ment are consistent with our results. In contrast, as we 
have already pointed out that our treatment is a non- 
perturbative one. Since any non-perturbative treatment 
is always welcome in field theory we thought it would be 
interesting to present these results and arguments. Be- 
sides we have been able to show the renormalisation of 
the free energy. These regularisation schemes help one 
to understand various interesting renormalization of the 
quantity involved (in this case charge) and a proper jus- 
tification of the scaling functions discussed in the paper. 

In summary, we have shown exactly, that charged 
scalar fields in (2+1) and (3+1) dimension at all tem- 
peratures are diamagnetic. 
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APPENDIX A: COMMENT ON AN ANOTHER 
REGULARISATION SCHEME 

In this appendix we would like to discuss about another 
rcgularisation scheme which at first sight one might think 
of as a correct regularisation scheme. However as we 
show below this method is not a correct scheme. This 
scheme can be thought of as relating L and A by the 
criterion that the maximum energy of the two systems 
should be matched. The matching of the energy cutoff 
gives the relation between the momentum cut off A and 
the Landau level L as 



E L = i/m 2 + (2L + l)eB = \J k 2 + m 2 (Al) 

This immediately suggests that 

A 2 = (2Z + l)e£? (A2) 

Note that this relation is different from the mode match- 
ing relation (p0|) we used previously for the charged scalar 
field theory. However, it turns out that this renormalisa- 
tion procedure gives a divergent answer for the free en- 
ergy difference A/. We support this statement through 
both numerical as well as analytical calculations. 

Below we present numerically a comparison between 
the two schemes. The values of Af A (scheme due to 
anonymous) and A/ J (used in the paper) are tabulated 
for B = 2, m = .1 and for various values of the cutoff L. 



Different values of cutoffs(L) Af J 



L=10 
L=100 
L=1000 
L=10000 
L=50000 
L=80000 
L=100000 



6.774 
20.308 
63.507 
200.248 
447.457 
565.929 
632.698 



0.2157 
0.2325 
0.2382 
0.2400 
0.2404 
0.2405 
0.241 



From this table it is easy to see that our method gives 
a finite convergent, cutoff independent answers while the 
other method gives cutoff dependent, divergent answers 
for the free energy. In fact one can show analytically in 
this case that A/ diverges with the cutoff as yL. This 
unphysical and meaningless answer is due to the com- 
parision of systems with different number of degrees of 
freedom. Note that in our analysis A/ converges as -7=. 
Hence, the scheme due to matching of highest energy is 
therefore clearly unphysical. However, as already dis- 
cussed our scheme can be justified physically from the 
adiabtically turning on the magnetic field. Then the en- 
ergy of each mode is affected by the magnetic field. Com- 
paring systems with the same number of modes is the log- 
ically and physically correct procedure. What the other 
method does is that it compares the systems with dif- 
ferent degrees of freedom. Hence, regularisation scheme 
based on energy cutoff does not give a renormalisable 
answer. Also notice that phase space is invariant in any 



Lorentz frame while the energy is not. We present be- 
low analytical work to show that in our case Sf varies 
with ^7= and as L — * 00, it does converge while the 

other case it diverges as vZ in the same limit as indi- 
cated above. Since the ultimate aim is to take the cutoff 
to infinity let us calculate the behaviour of the leading 
term in free enrgy difference as a function of cutoff. The 
difference between the free energy (using the two regu- 
larisation schemes) is given by 



5f L {B,m) 



1 r 



(m 2 + (2L+ \)eBf /2 
(m 2 + 2(L + l)eBf /2 



(A3) 



Now defining a dimensionless small parameter rj = 
(2L+i)eB (which is reasonable in large cutoff L and in 
strong magnetic field) it is easy to write the above differ- 
ence as 



1 



6f L (ri)~VL + Oi-j=) 



(A4) 



This is consistent with the numerical results presented 
in the above table. Another way of qualitatively under- 
standing the above divergence is the following . In a 
strong magnetic field limit all the higher Landau levels 
are far away and one is interested only in the lowest Lan- 
dau level. Now in this situation if one countes the energy 
states and matches with the free case, then there will be 
leading order diveregnce with the cutoff due to mismatch 
of number of degrees of freedom. Thus we notice that 
though the regularisation scheme provides a finite posi- 
tive free energy differenece, it fails to give a renormalised 
difference. Hence this scheme should not be considered 
as a correct and physical regularisation prescription. 



APPENDIX B: UNIQUENESS OF THE 
REGULARISATION SCHEME USED 

In this appendix we want to show the uniqueness of the 
regularisation method used in this paper. In particular, 
we would like to argue that the mode matching relation 
( pp| ) can be obtained from the condition that the differ- 
ence between the two free energies is independent of the 
cutoff and tends to some finite value. Here too we assume 
that the difference is finite for all values of B > and 
m > 0. As we have noted before for charged scalar fields 
the difference between the free energies can be written as 



A/ = eB V™ 2 + (2/ + l)eS - f pdp Vp 2 
1=0 J ° 



(Bl) 



We want to find the relation between A and L for which 
the difference tends to a finite value independent of the 
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cutoffs. In otherwords as L — > oo or A — > oo the differ- 
ence between the free energy becomes unique and finite. 
Though in principle A/ depends on both A and L sep- 
arately, we will notice that because of a unique relation 
between them we will get a finite value independent of 
both of them. Hence the free energy difference can be 
written as 

A/ = A/ + O(-^t), (B2) 

where the first term is the universal component which is 
independent of the cutoff (its value of course depends on 
eB and m) and the second term is the correction to it. 
The function g(L) is assumed to be a smooth function 
of L and in the present analysis we are not interested in 
the actual form of it. We assume the cutoff to be high 
enough so that we can replace the sum by an integral. 
Now if we demand that A/o should be independent of 
the cutoff then we get an equation of the form 

dA 



k—^m 2 +k 2 -► eB\/m 2 + 2LeB (B3) 
O Li 

We have used the symbol — > instead of an equal sign to 
signify "asymptotically approaching" rather than equal- 
ity between the two relations. In fact they are strictly 
equal only if the cutoff is taken to infinity. But, because 
of finite value of L or A however large, we do not expect 
strict equality in the above relation. From this equation 
in the first order approximation it is quite evident that 
there should be a relation between A and L such that 

This immediately gives a relation between the two cutoff 
within an undetermined constant C as 

A 2 = 2eBL + 2C (B5) 

It is important to notice that this relation has been ob- 
tained in the limit of large enough cutoff. In fact as we 
have seen in the previous subsection the relation is valid 
for small as well as large value L. However, the physical 
energy difference approaches to finite value as -4= as we 
vary L. In the next we will try to fix this constant from 
the fact that the difference goes to a finite value. We 
subst itut e this value in the free energy difference equa- 



tion (Bl) and plot the difference as a function of various 
values of L for fixed values of B and m. Then it turns 
out that the only allowed value of C for the free energy 
to be finite and independent of the cutoff is eB. Thus we 
find that the exact relation between the two cutoff is 

A 2 = 2eBL + 2eB = 2eB(L + 1) (B6) 

It is interesting to note that no other values of C are 
allowed for the free energy difference to be finite and 
independent of cutoff. Thus we have justified the mode 



matching relation (20). In fact, this is the only method 
by which one can get a unique finite cutoff-independent 
free energy difference. 
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